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I. INTRODUCTION 



Entanglement as the central feature of quantum mechan- 
(-h its distinguishes a quantum system from its classical coun- 
cil ^ .erpart. As an important physical resource, it has many ap- 
plications in quantum information theory. Among the well 
£3 known applications of entanglement are supcrdense cod- 
C^ing pj, quantum state teleportation@, Efforts to quan- 
^ tifv this resource are often termed entanglement theory 
, Q uantum entanglement also has many different applica- 
tions in the emerging technologies of quantum computing 
^ — I and quantum cryptography 0, Q , and has been used to re- 
^alize quantum teleportation experimentally Q- 
JO Quantum entanglement has attracted a lot of attention in 
recent years in various kinds of quantum optical systems @- 

■<^j- \ Several methods to quantify entanglement have been pro- 
posed. For pure states, the partial entropy of the density 
matrix can provide a good measure of entanglement. Infor- 
mation entropies are also used to quantify the entanglement 
;> in quantum information [T3 | . In this regard the von Neu- 
mann Entropy (NE) (la |. Linear Entropy (LE) and Shan- 
non information Entropy (SE) have been frequently used 
J^in treating entanglement in the quantum systems. It is 
worth mentioning that the SE involves only the diagonal 
elements of the density matrix and in some cases gives in- 
formation similar to that obtained from the NE and LE. On 
the other hand, there is an additional entropy, namely, the 
Field Wehrl Entropy (FWE) This measure has been 

successfully applied in description of different properties of 
the quantum optical fields such as phase-space uncertainty 
decoherence [HH3] etc. 



The FWE is more sensitive in distinguishing states than 
the NE since FWE is a state dependent (2l[ . The concept 
of the Wehrl Phase Distribution (WPD) has been devel- 
oped and shown that it serves as a measure of both noise 
(phase-space uncertainty) and phase randomization [2lj |. 
Furthermore, the FWE has been applied to the dynami- 
cal systems. In this respect the time evolution of the FWE 
for the Kerr- like medium has been discussed in [22] . For 
the Jaynes-Cumming model the FWE gives an information 
on the splitting of the Q-function in the course of the col- 
lapse region of the atomic inversion as well as on the atomic 
inversion itself [U, ES Hi- 

In the current contribution we study the evolution behav- 
ior of entanglements in a semiconductor cavity QED con- 
taining a quantum well coupled to the environment by 
the FWE, generalized concurrence vector and WPD. We 
also explore the situation in which the entanglement de- 
cays to zero abruptly. Recently, Yu and Eberly [25l— I20I] 
showed that entanglement loss occurs in a finite time un- 
der the action of pure vacuum noise in a bipartite state of 
qubits. They found that, even though it takes infinite time 
to complete decohcrence locally, the global entanglement 
may be lost in finite time. This phenomenon of sudden loss 
of entanglement has been named as "entanglement sudden 
death" (ESD). . Opposite to the currently extensively dis- 
cussed ESD, Entanglement Sudden Birth (ESB)P,iJ is 
the creation of entanglement where the initially unentan- 
glcd qubits can be entangled after a finite evolution time. 
These phenomena have recently received a lot of attention 
in cavity-QED and s pin chain|32l. [33j . and have been ob- 
served Exprimentally [H, [35| . 

The paper is organized as follows; Section 2 displays the 
physical system and its model Hamiltonian. Section 3 de- 
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votes to evolution equations of the system by the quantum 
trajectory approaching. Section 4 discusses the entangle- 
ment due to Wehrl entropy, generalized concurrence and 
the Wehrl phase distribution and section 5 supplies a con- 
clusion and outlooks. 



II. MODEL 

The considered system is a quantum well confined in a 
semiconductor microcavity. The semiconductor microcav- 
ity is made of a set of Bragg mirrors with specific separation 
taken to be of the order of the wavelength A. In the sys- 
tem under consideration, we restricted our discussion to the 
interaction of electromagnetic field with two bands in the 
weak pumping regime. The electromagnetic field can make 
an electron transition from valance to conduction band. 
This transition simultaneously creates a single hole in the 
valance band which leads to generation of exciton in the 
system. One can use an effective Hamiltonian without spin 
effects for describing the exciton-photon coupling in the 
cavity as [36U431 ]: 

H = hiOpO 1 a + hiuj^b + ihg 1 (a)b — b^a) 

+ ha'tftfbb + ih{e'e wt <J - h.c) + H r , (1) 

where u> p and uj e are the frequencies of the photonic and 
excitonic modes of the cavity respectively. The bosonic 
operators a and b are respectively describing the pho- 
tonic and excitonic annihilation operators and verifying 
[a,a^] = 1; [b,b^] = 1. The first two terms of the Hamil- 
tonian describe respectively the energies of photon and ex- 
citon. The third term corresponds to the photon-exciton 
coupling with a constant of coupling g' . The forth term 
describes the nonlinear exciton-exciton scattering due to 
coulomb interaction. Where a' is the strength of the inter- 
action between excitons 0, (45| . The fifth term represents 
the interaction of external driving laser field with the cav- 
ity, with e' and ui being respectively the amplitude and 
frequency of the driving field. Finally, the last term de- 
scribes the relaxation part of the main exciton and photon 
modes. We restrict our work to the resonant case where the 
pumping laser, the cavity and the exciton are in resonance 
{lo = ujp = uj e ). We have neglected also the photon-exciton 
saturations effects in Eq. (TT]). It is shown that these effects 
give rise to small corrections as compared to the nonlin- 
ear exciton-exciton scattering [13, 0, S3]- Furthermore, 
we assume that the thermal reservoir is at the T = and 



we neglect the nonlinear dissipations [481 ] . then the master 
equation can be written as |49l - l53| 

^ = -ia[bWbb, p] + g[{o)b - tfa),p] 

+ e[(a t -a),p]+Lp, (2) 

where t is a dimensionless time normalized to the round 
trip time r c in the cavity, and we normalize all constant 
parameters of the system to l/r c as: g = g'r c , e = e'r c , a = 
a'r c . Lp represents the dissipation term associated with 
H r and it describes the dissipation due to the excitonic 
spontaneous emission rate 7/2 and to the cavity dissipation 
rate k: 

Lp = n(2apa) — a) ap — pa' a) 

+ ~f/2(2bptf -tfbp- ptfb). (3) 

III. EVOLUTION EQUATIONS 

In the weak excitation regime - <?; 1, we can neglect the 
non-dia gon al terms 2apa^ and 2bpb^ in the master equa- 
tion® [HiL l55l]. The density matrix can then be factorized 
as a pure state [45|,[54|]-[53- We then obtain, the following 
compact and practical master equation: 

t = MHeffP-(H effP n (4) 

where the effective non-Hcrmitian Hamiltonian H e ff de- 
fined as 

H eff = ihg(a)b - 6+ a) + hatftfbb 

+ ihe(a < — a) — ihna'a — ih—b^b. (5) 

in which the time dependent density matrix p = 
\ip(t))(ip(t)\ is a possible solution of equation Also \ip(t)) 
satisfies the following equation: 

lh mm = Heffm)Y (6) 

The essential effect of the pump field is to increase the 
excitation quanta number in the cavity which allows us to 
neglect the term hea in the expression of the effective non- 
Hermitian Hamiltonian Eq.© [54l - l5(| . 
We can expand \ip(t)) into a superposition of tensor product 
of pure excitonic and photonic states [45| . [5^ - l56j : 

|V (i)) = Ax)|00) + A w \10) + AnlOl) + Aulll) 
+A 20 |20) + A 02 |02) + A 30 |30) + A 03 |03) 
+A 21 |21)+A 12 |12), (7) 
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where \ij) = \i) (^) \j), is the state with i photons and j 
excitons in the cavity. We then obtain the following differ- 
ential equations for the amplitudes Aij(t) 



Shannon entropy Sh depends on the diagonal elements so 
that it does not contain any information about the phase 
and can be expressed as 



dA 00 

A 

dAu 

As 

dA 3 Q 

d%i 
dt 



-eA w , 

-eAu - gA 10 - 2^01, 



e (A 0Q - V2A 2a ) + gA 01 - kA w , 

V2g (A 02 -A 20 )-(k+ |) A n + sA 01 - ey/2Aai, 

V2gA n ~ 2kA w + V2eA 10 - V3sA 30 , 

-V2gAn - 2iaA 02 - 7^02 - eAu, 

~gV3A 12 - (^ + 6ia)A 03 (i), 

eV3A 20 + gV3A 21 - 3kA 30 , 

eA 02 + g {V3A 03 - 2A 21 ) - (k + 7 + 2ia) A 12> 



eVZAn + g (2A 12 - V3A 30 ) 



2k 



A 21 , 



(8) 



We assume that, at time t = the vector state \ip(t)) is 
in vacuum state, \tp(t = 0)) = 1 00) : 



A ii {t = 0) = 0. 



(9) 



For pure state, the density operator can be written in term 
of the wavefunction \ip(t)) as p p h.exc = IVKOK^MI- The 
reduced density matrices of photon-exciton system can be 
written as 



P P h = tr exc (\iP(t))(ij(t)\),p exc = tr ph (\^(t)){^(t)\). 

(10) 

the above equation will be used in the next sections exten- 
sively to calculate the FWE, concurrence and WPD. 



IV. ENTAGLEMENT DYNAMICS OF THREE 
EXCITATIONS REGIME 

A. Wehrl entropy 

In this section, we investigate the field Wehrl entropy 
for the system under consideration. Actually, the Wehrl 
entropy is better than the Shannon entropy and von Neu- 
mann entropy for certain states. More illustratively the 



S H {t) = -Jj>(n,t) Hp(n,t)), 



(11) 



where p(n, t) = (n\p (t) \ n) is the photon number distribu- 
tion. On the other hand, the von Neumann entropy defined 
as S]\r(t) = — Tr(p(t) ln(p(i)), can not be used in the mixed 
state case. 

To study of the Wehrl entropy of the photons in the case 
of three excitations regime one need to calculate the Husimi 
Qf function. Which is defined in terms of the diagonal 
elements of the density operator in the coherent state basis 
as 



} ph {p,®,t) = -\{P\i>{t)){' 

IT 



(12) 



where the coherent state representation \0) = b n {0) \n) 

n=0 

while the amplitude 



bnifi) = ^ J 8 ' V , 



(13) 



In order to compute the Husimi Q function of the photons 
we substitute the sate vector in the case of three photon 
excitation is given by Eq. {J7J) into Ea. (fl"2"l) which reads 



? P h03,e,t) 



qoW)A 00 + qx{p)A 1Q + q 2 (f3)A 20 + q 3 {/3)A 30 

2 

q W)A 01 +(&(P)Au +q 2 ((3)A 21 

2 

qo(P)A 02 + qi (f3)A 12 + q Q (f3)A 03 



(14) 



Now, we may calculate the Wehrl entropy. The concept 
of the classical-like Wehrl entropy (FWE) is a very infor- 
mative measure describing the time evolution of a quantum 
system. The Wehrl entropy, introduced as a classical en- 
tropy of a quantum state, can give additional insights into 
the dynamics of the system, as compared to other entropies. 
The Wehrl classical information entropy is defined as |16j 

/>27T />OG 

Sw(t) = - / Q p h (A 0, t) In Q ph ((3, 0, t) \/3\ d \/3\ dQ 
Jo Jo 

(15) 
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(k,g,e)=(0.00002,0.1,0.2) a10 -a ( T ,g,e)=(0,0.1,0.2) do" 8 




FIG. 1. The time evolution of the Wehrl entropy Swif) f° r 
a = I0~ 8 , (fe,g,e)=(0.00002,0.1,0.2) and with different values 
of 7. 

We point out that the state vector coefficients and Q- 
function both are normalized at time steps as follows 

£MM 2 = 1 (16) 

/ / Q ph (P,e,t)\P\d\0\dG = l (17) 
Jo Jo 

To explore the influence of decoherence on the dynami- 
cal behavior of the Wehrl entropy, we have plotted the time 
evolution of the photon Wehrl entropy Sw (t) as a function 
of time t for different values of the coupling constant 7 and 
the cavity dissipation rate k in Figures. (1) and (2) 
Figure. (1) shows the influences of excitonic spontaneous 
emission rate 7 on the Wehrl entropy (WE). By increas- 
ing 7 the (WE) decreases. Furthermore, similar effect for 
the cavity dissipation rate k can be observed in the Fig. (2). 
It is worth to note that for large values of k the (WE) 
decreases abruptly much faster than Fig.(l). The increas- 
ing of k or 7 enhances the decoherence in the system and 
consequently causes the destruction of entanglement in the 
system. To have further insight, we plot in Fig. (3) and 
Fig. (4) the Wehrl entropy for different values of the cou- 
pling constant g and the amplitude of the driving field 
e, respectively. By increasing the coupling constant g the 
frequency oscillation of the Wehrl entropy increases . This 
effect is also observed in the autocorrelation function [59| 
and in two photon excitations [fjoj . 



FIG. 2. The time evolution of the Wehrl entropy Sw(t) for 
a — 10~ 8 , (7, g, e)=(0, 0.1, 0.2) and with different values of k. 
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FIG. 3. The time evolution of the Wehrl entropy Sw(t) for 
a = 10" 8 , (k, 7, e)=(0. 0002,0,0. 2) and with different values of 
9- 



B. Generalized concurrence 

To study of entanglement for pure states usually the par- 
tial entropy of the density matrix is a good measure of en- 
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FIG. 5. The time evolution of the Concurrence C(t) for a = 
10" 8 , (fc,ff,e)=(0.0002,0.1,0.2) and with different values of 7. 

tanglcment which reads 

E(tp) = -tr(p Lnp) = - £ (A 4 LnXi), / lg \ 

i 

where p is the reduced density matrix, A^ is the i th eigen- 
value of p. In the case of a two qubit mixed state p, the 
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FIG. 7. The time evolution of the Concurrence C(t) for a — 
10~ 8 , (k, 7,e) = (0.0002,0,0.01) and with different values of g. 

concurrence of Wootters can be used as a measure of en- 
tanglement which is given bvfolj 

C{p) = max(0, Ai - A 2 - A 3 - A 4 ), (19) 

in which the Ai are the square roots of eigenvalues in de- 
creasing order of ^fpp^fp with p = {a y (^) a y )p*{a y &)o y ). 
Recently, some extensions have proposed for definition 



G 



of concurrence in the case of an arbitrary bipartite pure 
state lib) = 



(b) T-0-82 



J2 £ a i] \ e i ® e i) as |63,|6i 



i=l j=2 



/ATI iV2 



v v »<j fe <* 



(20) 



where N = min(Ni, N 2 ). 

Here, we deal with a pure state £ C 4 (^)C 4 so that, 
N = Ni = N 2 = 4. To study the time evolution of the 
concurrence in the case of three excitations we substitute 
the state vector ((?]) into the Eg. (|2T))) . thus we obtain 



(21) 



where T(t) is given by 

T(t) = \A 00 A 21 - A Q1 A 20 \ 
+ |A 00 Aii-^oiAio| 
+ \A 01 A 12 - A 11 A 02 \ 
+ \A 00 A 12 - A 12 A W \ 
+ \A 10 A 21 - AuA 20 \ 
+ \A 12 A 20 \ 



AosAwl 
A 03 A 12 \ 

^20^03! 

A 03 A 21 \ 
^4o2^3o| 

A 2 iA 12 \ 
A 12 A 30 \ 



A n A 03 \ 
Ao 2 A 2 o\ 
A 02 A 21 1 
An -430 1 

-4o3^3o| 



A U A 
A 21 A 



(22) 



We plot the time dependent concurrence vector as a func- 
tion of time t for three values of 7 in Fig. (5). As it seen 
any increasing of 7 leads to decreasing of entanglement sim- 
ilar to the Wehrl entropy. Furthermore, an interesting cases 
are observed in the Fig. (6)- (7). These figures show that the 
concurrence is periodic in the domain of time. Moreover, 
unlike the large values of e, figure (6) shows that entan- 
glement can fall abruptly to zero (the two lower curves in 
the figure) for small values of e(e — 0.025 and e = 0.02), 
and remains zero for a period of time before entanglement 
recovers. The abrupt disappearance of entanglement that 
persists for a period of time is referred to as sudden death 
of entanglement (ESD)[25|, [2(| |fjH |65[ and also the fast 
appearance of entanglement after a while is called sudden 
birth of entanglement(ESB) [H El . The length of the time 
interval for the zero entanglement is dependent on the val- 
ues of e. The smaller values of e, the longer the state will 
stay in the disentangled separable state. Furthermore we 
show that the (ESD) and (ESB) can be affected strongly 
by the coupling constant g. As it is seen from figure. (7) 
the (ESD) and (ESB) can be enhanced by increasing of 
coupling constant g. We point out that, in the Figures. (5)- 
(7) we assume the zero value of the excitonic spontaneous 
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FIG. 8. The time evolution of the Wehrl PD Se(t) for a = 1CT 
(k, g, e)=(0. 0002,0,0. 2) and with different values of 7. 
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FIG. 9. The time evolution of the Wehrl PD S&(t) for a = 10" 
(e, g, 7) = (0.2,0,0) and with different values of k. 



emission rate 7 = 0, thus one important reason for the 
(ESD) is the interaction of system with its surrounding. 
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C. Wehrl Phase distribution 

The Wehrl phase distribution (Wehrl PD), defined to be 
the phase density of the Wehrl entropy [l7|, HH], i.e., 

Se(t) = -[ Q ph (f3,Q,t)\nQ ph (p,e,t)\l3\d\f3\, (23) 

where 6 = arg (/?) and Q p h (/3, 6, t) is given in Eq. (fT2"j) . 

Based on Eq. (|2"5)) . we present some interesting results 
for the effects of excitonic spontaneous emission and the 
dissipative rate of the cavity on the entanglement behavior 
in the point of view of Wehrl PD. It is observed that when 
7 = (see fig. 8(a)) Se(t) oscillates between maximum and 
minimum peaks which is an indication of ESB and ESD. 
For 7^0 the situation is completely different, the excitonic 
spontaneous emission destroys the entanglement (see fig. 8). 

Now, we would like to answer the question: How Se(£), 
is influenced by the cavity dissipation? For this purpose, 
we take two different values of k in fig. 9. For small values 



of k S@(t) oscillates but when k increases Se(*) decreases 
quickly without oscillation(see figure 9(b)). This shows a 
one-to-one correspondence between the behavior of S&(t) 
and the Wehrl entropy or concurrence which opens the door 
for using S@(t) as an entanglement measure. 



V. CONCLUSION 

In this paper we have studied the dynamical behavior 
of the quantum entanglement for a semiconductor micro- 
cavity containing a quantum well. The system is pumped 
with weak laser amplitude. We studied the time evolution 
of entanglement between the photon-exciton by the field 
Wehrl entropy, generalized concurrence and Wehrl phase 
distribution. Our results show that the new features such 
as entanglement sudden death and entanglement sudden 
birth can be reported for specific values of the cavity dissi- 
pation rate and the excitonic spontaneous emission rate. 
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